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Abstract 

We analyze the shapes of roughness distributions of discrete models in the Kardar, Parisi and 
Zhang (KPZ) and in the Villain, Lai and Das Sarma (VLDS) classes of interface growth, in one 
and two dimensions. Three KPZ models m. d = 2 confirm the expected scaling of the distribution 
and show a stretched exponential tail approximately as exp (— x*^'^), with a significant asymmetry 
near the maximum. Conserved restricted solid-on-solid models belonging to the VLDS class were 
simulated in d = 1 and d = 2. The tail in d = 1 has the form exp {—x^) and, in d = 2, has a simple 
exponential decay, but is quantitatively different from the distribution of the linear fourth-order 
(Mullins-Herring) theory. It is not possible to fit any of the above distributions to those of 1//°^ 
noise interfaces, in contrast with recently studied models with depinning transitions. 

PACS numbers: PACS numbers: 05.40.-a, 05.50.+q, 68.55.-a, 81.15.Aa 
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I. INTRODUCTION 



The simplest quantitative characteristic of an interface is its width, or roughness, w, 
defined as the rms fluctuation of the height around its average position. The scaling of the 
average width of different realizations with the length L and the time t is useful for analyzing 
interfaces formed in several processes . Another potentially useful characteristic is the 
full roughness distribution at the steady state, particularly in cases where reliable estimates 
of scaling exponents are not available. 

The distributions for linear growth models, such as the Edwards- Wilkinson (EW) and 
the Mullins-Herring (MH) equations, were obtained exactly in d = 1 and can be computed 
with high degree of accuracy in d = 2 P, 4, |^. The generating functions for Gaussian 
interfaces with power spectrum of the type 1//" were computed by the same methods 
However, many real interfaces are described by nonlinear equations at a coarse-grained 
level, and their scaling exponents and roughness distributions are exactly known only in 
some particular cases. Consequently, their accurate numerical calculation may provide a 
basis for comparison with results from experiments or model systems. In this Brief Report, 
we will study numerically the distributions of models governed by second and fourth order 
nonlinear growth equations, in d = 1 and d = 2. 

n 

The first model is that of Kardar, Parisi and Zhang (KPZ) jSj], in which the local height 
h at position x and time t evolves as 

— = y,y^h + \2{Vhf + r^{x,t). (1) 

Here, V2 and A2 are constants and 77 is a Gaussian noise with zero mean and variance 
(?7 (x, t) ri{x' ^ t')) = D6'^{x — x')6 (t — t'), where d is the dimension of the substrate. In the 
case A2 = 0, we obtain the EW equation. 

Another important nonlinear model is that of Villain and Lai and Das Sarma the 
so-called VLDS equation (or nonlinear MBE equation) 

— = u.V'h + KV\Vhf + r^(f , t), (2) 



where 1/4 and A4 are constants. The linear version of Eq. (0) (A4 = 0) is the MH equation 
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3. 



n 



Previous studies j3|, |5|] of the roughness distributions of discrete models in these classes 
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suggested that they have a universal shape, with the form 

Pl (w,) = -i-<l> (^] , (3) 

where Pl (^2) dw2 is the probabihty that the width of a given configuration hes in the range 
[w2, W2 + dw2\. Alternatively, one may use the scaling form 

= (4) 

(J \ a 



where (t = \J {w2^) — (^2)^ is the rms deviation of the squared width. 

For KPZ m d = 1, these scaling relations are analytically predicted and confirmed by 
numerical data j^. In d = 2, simulations of discrete models showed agreement with the 

was not analyzed in detail. 

Here, we will calculate accurate roughness distributions for three discrete models in d = 2: 
the etchin g ra odel of Mello et al Q], the restricted solid-on-solid (RSOS) model of Kim and 
Kosterlitz and a generalized RSOS model with maximum neighboring heights difference 
Ai/^,, = 200, hereafter called RS0S2. Their analysis is presented in Sec. II. 

Roughness distributions of the VLDS class were previously obtained by Racz and Plis- 
chke p|, who simulated a model with Arrhenius dynamics m. d = 2. In order to analyze 
the shapes of the VLDS distributions, both in = 1 and d = 2, we will simulate conserved 
solid-on-solid (CRSOS) models which show moderate finite-size effects when compared to 



other discrete VLDS models. Their original versions we re p roposed by Kim et al [17[ and 
were rigorously proved to belong to the VLDS class [l^IlMl- Their analysis is presented in 
Sec. III. 

In Sec. IV we summarize our results and conclusions. 



II. ROUGHNESS DISTRIBUTIONS OF KPZ MODELS IN 2 + 1 DIMENSIONS 

In Fig. 1 we show the distributions of KPZ models scaled according to Eq. with 
X = {w2 — {102)) jo. Results for other system sizes confirm the expected scaling of Fig. 1. 
The number of configurations in which W2 was measured was near 10^ for each model and 
periodic boundary conditions were adopted. Eq. Q and (jH) are equivalent, but the latter 
was used because it provides much better data collapse for the KPZ models. 
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Quantitative information on the shape of the KPZ distribution can be obtained 
from the dimensionless ratios of the moments of the roughness distribution Wn (L) = 
Io° ("^2 — {u!2))"'Pl The ratios obtained with better accuracy involve the lowest-order 
moments: the skewness 5* = ^3/2 and the kurtosis Q = — 3. Although these quantities 
do not characterize the whole distribution, they may be useful for the comparison with data 
from experiments or other model systems if data collapse methods are unreliable. 

The finite-size estimates of 5* and Q for the distributions of the RSOS and the etching 
models are shown in Figs. 2a and 2b. Their asymptotic estimates are S = 1.70 ± 0.02 and 
Q = 5.4 ± 0.3. These large values show that the roughness distribution strongly deviate 
from a Gaussian near the maximum. For x < xm, where xm is the abscissa of the maximum 
of the distribution in Fig. 1, the curve has an approximately Gaussian shape centered at 
Xm- However, for a; > xm, it significantly deviates from a Gaussian, showing a much slower 
decay. 

The tail of the log-linear plot in Fig. 1 is not a straight line, contrary to what is observed 
in the 1 + 1-dimensional KPZ model. This is certainly related to the non-Gaussian behavior 
of the KPZ interface .2.1 d„„e,.,o. Q Q. It . .ea.o„ab,e to a.u.e that the .aUng 

function decays as (x) ~ x^exp {—Ax'^). Thus, for fixed /5, the estimates of the exponent 
7 are given by 

In 



7(x) 



ln(x-'3*{x')) 



ln((x-A)~'^<^(x-A)) 



(5) 



ln[a;/(2;- A)] 
with constant A. 

In Fig. 3 we show the effective exponents obtained from the data of the RSOS and the 
etching models in L = 256 as a function of 1/x^, for three values of (3. The assumption of 
pure exponential decay (/5 = 0) leads to an asymptotic estimate 7 ~ 0.8. With P = —0.2, we 
obtain the smallest fluctuation of 7 (x) in the range of the variable x analyzed here. Results 
for values of (3 near —0.2 give 7 = 0.8 ±0.1. Although j3 = —1 suggests a simple exponential 
decay (7 = 1) for the etching model, it shows discrepancies with the estimates for the RSOS 
model, reinforcing the proposal of the stretched exponential form. 

One question that could be raised here is the possibility of fitting the KPZ distributions 
by those of 1//" signals However, the plots of the distributions for several values of 
a show that this procedure does not provide good fits of our data. Indeed, this can be 
quantitatively explained by the above values of the skewness and the kurtosis, which are 
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much larger than the upper bounds for the distributions of 1//" signals [3]: S = -\/2 and 
Q = 3, obtained in the limit a oo. Moreover, the simple exponential decay of the 1//" 
distributions for large x disagrees with the exp (— x^ *^) decay of the KPZ models. 

III. ROUGHNESS DISTRIBUTIONS OF VLDS MODELS 

The CRSOS models considered here are extensions of the original ones in which the 
incident particle executes a random walk among neighboring columns until finding a position 



in which its aggregation is allowed j2^. In CRSOSl and CRS0S2 models, the conditions 
^Hmax = 1 and AHmax = 2 between neighboring columns are satisfied, respectively. Nearly 
10® configurations were used to measured W2 in each case. 

Here, we will show distributions scaled according to Eq. 0, as done in most previous 
works. This is also helpful for a comparison with known results for the MH theory. 

In Fig. 4a we plot log [(u;2)^l(w2)] versus y = ^ for the CRSOSl and the CRS0S2 
models in d = 1, with L = 256. The scaling function for the MH theory, given in Ref. Q], is 
also shown. Its simple exponential decay contrasts with the more rapid decay of the VLDS 
models. 

The skewness and the kurtosis of the scaled VLDS distribution are obtained along the 
same lines of the KPZ values: 5* = 1.22±0.05 and Q = 1.6±0.2. They contrast with the much 
larger values S = 1.988 . . . and Q = 5.951 ... of the MH theory [J, Q]. An approximately 
exp (— y^) decay is suggested by a linear fit of the data for the CRSOSl model in the range 
2 < y < 5, as shown in Fig. 4b. 

Our results in d = 2 are shown in Fig. 4c, where scaled distributions of the CRSOSl 
model (L = 64) and the MH theory are compared. The latter was obtained from the 
generating function given in Ref. Q and the calculation of the residues of more than 20 
poles closest to the origin. Both distributions have a simple exponential decay for large y, 
which is justified in the VLDS case by the linear fit shown in Fig. 5 (2 < x < 4). However, 
the VLDS distribution is much narrower than that of the MH theory: the maximum of the 
VLDS scaling function is $ ~ 1.3, while the maximum of the MH one is $ ~ 1 (see, e. g. 
Ref. j,5|), and the respective decays are, approximately, exp(— 1.8y) and exp (—1.3?/). For 
the VLDS distribution, the skewness and kurtosis are 5* ~ 1.1 and Q ^ 1.8, while the values 
for the MH theory are S" ^ 1.30 and Q ^ 2.63. 
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ed distributions of the CRSOSl model, in d = 1 and d = 2, 
3], with some value of a. In = 1, no good fit was already 



We also tried to fit the sea. 
with l//"-noise distributions 
expected due to the form of the decay of the scaling function discussed above. In d = 2, 
there is no simultaneous agreement of the values of S and Q of those distributions and our 
estimates, for any a, indicating that this type of fit is also impossible. In Ref. P], a fit with 
a = 3 was suggested for the distribution of a growth model with Arrhenius dynamics, but 
this was probably a consequence of the lower accuracy of the data for that model. 

IV. CONCLUSION 

Roughness distributions were calculated for several discrete models in the KPZ and the 
VLDS class, in one- and two-dimensional substrates, confirming the scaling relations Q and 
(0)). The scaling functions are different from their counterparts of the linear theories, EW 
(second order) and MH (fourth order). In the 2 + 1-dimensional KPZ class, large values of 
the skewness 5* and kurtosis Q of the scaling functions are related to the high asymmetry 
near the maximum. For large x, the scaling functions decay approximately as exp {—x^'^), 
in contrast to the simple exponential decay of the linear theories and of 1 + 1-dimensional 
KPZ. In the 1 -|- 1-dimensional VLDS class, S and Q are smaller than those of the MH theory 
and a decay approximately as exp (— x^) is obtained. For the 2 + 1-dimensional VLDS class, 
the simple exponential decay of the scaled distribution is quantitatively different from that 
of the MH theory. 

Another interesting point of our study is that none of the above distributions are fitted by 
those of 1//" noise, for any a. It contrasts with the distributions for model and experimental 
interfaces in the depinning thresholds analyzed in recent works 0, 13; 2^ , which indicates 
strong effects of the non-Gaussian nature of the KPZ (in d = 2) and VLDS interfaces. On 
the other hand, it suggests the study of other models with depinning transitions, particularly 
those with relatively simple dynamics whose growth phases are in the KPZ class [26 1, for a 
comparison with 1//" noise distributions. 
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FIG. 1: Scaled roughness distributions (according to Eq. 0J) at the steady states of the models: 
RSOS in L = 256 (solid line), etching in L = 256 (squares) and RS0S2 in L = 128 (triangles). 
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FIG. 2: (a) Skewness and (b) kurtosis of the interface width distributions at the steady states of 
the RSOS (squares) and etching (triangle) models as a function of inverse lattice length. 
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FIG. 3: Effective exponents 7 (x) of the exponential decay of the scahng functions versus for 
three values of the exponent f3 and the models RSOS (squares) and etching (crosses). The dotted 
(dashed) curve is a linear fit of the data of the RSOS (etching) model. 
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FIG. 4: (a) Scaled roughness distributions (according to Eq. ^ at the steady states of the CRSOSl 
(squares) and CRS0S2 (crosses) models in d = 1, with L = 256. The solid curve is the scaled 
distribution of the MH theory in d = 1. (b) The same distributions for CRSOSl and CRS0S2 
models as a function of and a linear fit (dashed line) of the data for the CRSOSl model, (c) 
Distributions of the CRSOSl model m d = 2,l#ith L = 64, and a linear fit of the data (dashed 
line) in the range 2 < y < 4. The solid curve is the distribution of the MH theory in c? = 2. 



